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Since the main experimentally testable prediction of grand unified theories is the instability 
of the proton, precise determination of the proton lifetime for each particular model is desir¬ 
able. Unfortunately, the corresponding computation usually involves theoretical uncertainties 
coming e.g. from ignorance of the mass spectrum or from the Planck-suppressed higher¬ 
dimensional operators, which may result in errors in the proton lifetime estimates stretching 
up to several orders of magnitude. On the other hand, we present a model based on SO(IO) 
gauge group which is subsequently broken by a scalar adjoint representation, where the leading 
Planck-suppressed operator is absent, hence the two-loop precision may be achieved. 


1 Introduction 

Baryon number violation as a smoking gun of the grand unified theories (GUTs) is searched for 
by number of experiments. Let us concentrate for simplicity on the p —> e + 7r° decay which has 
a clean signature in the water Cherenkov detectors and which is usually predicted by non-SUSY 
GUT models. For this channel the recent bound on the proton lifetime is set by the Super- 
Kamiokande experiment 1 to r p > 1.4 x 10 34 y, whereas the future Hyper-Kamiokande experiment 
assumes to reach the bound of t p > 10 35 y after 10 years of data taking 1 . If this progress on 
experimental side should distinguish between different GUT models, the uncertainties in the 
theorists’ computations of proton lifetime can not exceed one order of magnitude and a next- 
to-leading order computation is necessary. 

The aim of this text is twofold. First, we would like to explain in Section 2 what are the 
ingredients for such a higher order computation of the proton lifetime. Second, we would like 
to present a particular GUT model based on SO(IO) gauge group which was analyzed at next- 
to-leading order and interesting correlations between the proton lifetime and particle spectrum 
were found (see Section 3). 

2 Error estimates for proton lifetime computation in GUTs 

In non-SUSY GUTs the proton decay is predominantly mediated by vector bosons with large 
masses Mx and if the the contributions of the scalar fields are neglected, My is also the scale 
above that the /3-functions for the SM couplings coincide. For this reason, we assume in the 
first approximation that the scale Mq where the SM couplings intersect is equal to Mx, hence 
determines the proton lifetime r p . The errors in the position of Mq then lead to errors in r p , let 
us examine what are their sources. 

1. Finite order of perturbation theory. Defining the variable t = ^ log the renormaliza- 


tion group equation for the gauge couplings may be rewritten in the form 


d 

df 



— 



Cijk 

167T 2 


Ot-jOl-k + 


( 1 ) 


where cti = gf / 47T, a, b and c are 0(1) coefficients corresponding to one-, two- and three- 
loop contributions, and the dots correspond to higher order contributions. Consequently, 
one can estimate that if working at one-loop level, the error coming from neglecting the 
two-loop effects at the GUT scale M G ~ 10 16 GeV is of the order of 

Ac>! 1 “ 1 (Mg)i-Ioop ~ ~~ tz ) ~ 0,03 x ( 2 ) 

Similarly, the error of the two-loop calculation reads 

Aa l “ 1 (M G ) 2 _ioo P ~ (16^2^2 - tz ^> ~ °- 0002 x ^(1) • ( 3 ) 


2. Measurement of as at the EW scale. Taking into account the experimental value for the 
strong coupling as(Mz) = 0.1185 ± 0.0006, ler uncertainty in ad) 1 reads 

Aag 1 (Mc)ia ~ Aag 1 (M^)i cr = 0.04 (4) 

(strictly speaking, the uncertainty is reproduced at the GUT scale without any change 
only in case of one-loop running which is linear). Typically, the 0(1) factor in Eq. 2 
makes the one-loop error more significant than the error in Eq. 4, however, it is obvious 
that switching to the three-loop calculation is meaningless with today’s precision in as 
measurement. 


3. Threshold effects. In real models the masses of the heavy fields are not equal and so called 
matching- of the couplings at the GUT scale has to be performed. If a simple gauge group 
G is spontaneously broken into a direct product of subgroups Gi (with at most one abelian 
factor) at the scale )jl, then at the one-loop formula reads 5 
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(V, F and S denote the heavy vector bosons, fermions and scalars integrated out at the 
scale //). As a rule, all the heavy masses are products of the scalar couplings and the 
vacuum expectation value of the scalar field responsible for the breaking of the group G, 
and the matching scale g is chosen close to the barycenter of these masses. Consequently, 
the heavy masses M should lie close to /i and comparing Eqs. 5 and 2 one obtains 


A a i 1 (M G ) matc hing ~ Acq 1 (Af G )i_i oop . 
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~ log However, the threshold effects may be even more significant for 
itting of the heavy spectrum. Therefore, if one would like to work at two-loop 


precision level, the (one-loop) threshold effects have to be taken into account. 


4. 


Gravity induced operators. Since the unification scale is typically rather close to the (re¬ 
duced) Planck scale Mp\ = 2.43 x 10 18 , one has to take into account also the effective 
operators such as 1 
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Figure 1: On the l.h.s., the effect of A a. i 1 (Mg)i-ioo P on the determination of Mg is shown whereas the errors in 
the proton lifetime due to Aa“ 1 (MG)i-ioo P , Aag 1 (Mo)ia, and Aq“ 1 (Mg)pi are depicted on the r.h.s. 


where G pu is the GUT field strength, H is a scalar multiplet and k is an 0(1) constant. Let 
us suppose that the vacuum expectation value (VEV) (H) = Vq breaks the gauge group G 
to the SM. Then Vq ~ Mq and Eq. 6 contains the contributions to the SM kinetic terms 
e,;Tr (F^F^J with e* ~ 0( 10~ 2 ). In order to maintain the canonical normalization of 

the kinetic terms, one has to redefine A l p —> (1 + ej) 1//2 A p , g% —> (1 + e*) -1 / 2 ^. It is this 
redefined coupling which is measured at the EW scale, however, the unification condition 
holds for the original couplings: ac = (1 + ei)ai(Mc) for all i = 1,2,3. Naturally, this 
induces the uncertainty 

Acq 1 (A1g)v\ ~ 0(10 2 ) 

which is comparable to the one-loop error, let alone the fact that in case of large number 
of fields in the theory, the effective Planck scale may be lowered by as much as one order 
of magnitude 1 and then the “gravity smearing” effects are even more pronounced. 

To demonstrate the effect of the uncertainties in a” 1 on the proton lifetime prediction, let 
us consider a toy model where only the SM fields are present up to the unification scale Mq 
which is determined by the intersection of the SU (2 )l and strong couplings (this could be the 
case e.g. in the simplest flipped SU(5) models). Assuming Mx = Mq let us evaluate the proton 
lifetime using the naive formula 0 r p ~ 5 . If the uncertainties A ai(Mc)~ l are applied, the 

corresponding errors in Mq are found to be rather large due to the small angle between the 
two curves describing the running. On the l.h.s of Figure 1 we show an example of the effect 
of Aa^ 1 (Mg)! -loop on the determination of Mg (assuming the 0(1) factor in Eq. 2 to acquire 
values between —5 and 5), the corresponding errors in the proton lifetime are then shown in 
the first column of the plot on the r.h.s. The second column in the same plot corresponds to 
the proton lifetime computed using the two-loop /3-function with the error-bars depicting the 
lcr uncertainty in as(Mz )• Finally, the third column shows the error coming from the Planck- 
suppressed operators where we assumed k £ (0,0.5), > 0 and eg < 0 which leads to a shift to 

lower Mg and hence only one-sided error bar in r p . For Mg shifted to higher values, the gravity 
smearing effects get completely out of control and the two couplings may not even intersect. 

Finally, let us add that the 30-40% error in the lattice computation of the hadronic matrix 
elements 1 introduces also an uncertainty in r p , which, however, does not exceed one order of 
magnitude. On the other hand, the ignorance of the superpartner spectrum in case of the SUSY 
GUTs may lead to an error stretching up to several orders of magnitude. 

“We are aware of the fact that in the precise calculation the proper formula has to be used, however, we believe 
that this simplification is enough for estimating the errors which are mostly driven by the fact that r p oc Mq 















3 Minimal SO(IO) 


Fortunately, the gravity smearing effects are absent in some of the unification models such as 
the one based on SO(IO) gauge group broken by the VEV of a scalar adjoint representation 45// 
since Tr (G Mi ,45//G !Mi ') = 0, and the next-to-leading order analysis is possible. 

Although this model was abandoned for a long time due to the presence of tachyonic in¬ 
stabilities in the tree-level spectrum, it was found 1 that this problem is cured at the quantum 
level. Subsequently, the SO(IO) unification with scalar sector composed of 10-, 45- and 126- 
dimensional representations was studied in depth and it was observedthat the exact unification 
and also the seesaw scale a = (126) in the ballpark of 10 13 — 10 14 GeV may be achieved, if the 
scalar fields with SM quantum numbers (8, 2, +1/2) or (6, 3,1/3) are accidentally light. In both 
cases the scalar spectrum was computed, hence the two-loop analysis including the one-loop 
matching was performed"’’ and an interesting feature of the parameter space compatible with 
the unification and all the other constraints was found: Either the proton decay will be seen in 
the Hyper-Kamiokande experiment or only the light octet scenario is viable and this scalar is 
within the reach of LHC 5 . 

4 Conclusions 

Some of the obstacles to the precise determination of the proton lifetime in the GUTs were 
mentioned in this text. Namely, we show that the recent error in the measurement of the strong 
coupling does not allow the computation beyond the two-loop precision level. Moreover, the 
Planck-suppressed operators or ignorance of the mass spectrum may introduce an uncertainty 
comparable with the one-loop error. In this case even the two-loop precision can not be reached 
and using a toy model we show that the error in the proton lifetime estimate considerably 
exceeds one order of magnitude, the size of the improvement of the planned experiments with 
respect to the recent bounds. On the other hand, we presented an SO(IO) model where the 
Planck-suppressed operators are absent because of the group-theory structure, hence the two- 
loop precision may be achieved. 
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